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Abstract 

Solutions of a 1-D free-interface problem modeling solid combustion front propagat- 
ing in combustible mixture with periodically varying concentration of reactant exhibit 
classical phenomenon of mode locking. Numerical simulation shows a variety of locked 
periodic, quasi-periodic and chaotic solutions. 
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1 Introduction 

The current letter is intended to communicate new observations based on numerical solutions 
of the two-phase Stefan problem with kinetics. The simulations reveal some previously 
unknown features of this dynamically diverse system. Namely, it exhibits the phenomenon 
of frequency locking in response to a spatially-periodic perturbation of the medium. 

The free-boundary problem that is the subject of the paper arises naturally as a mathe- 
matical model of a variety of exothermic phase transition type processes, such as condensed- 
state combustion (also known as Self-propagating High-temperature Synthesis or SHS [313), 
solidification with undercooling ji], laser induced evaporation rapid crystallization in thin 
films [H^ etc. 

This study represents a natural extension of the numerical experiments described in pp, 
where it was demonstrated that, due to the competition between the heat release at the 
interface and the heat dissipation by the medium, the system generates a variety of complex 
thermokinetic oscillations. The dynamical patterns exhibited by the unperturbed system. 
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as the governing parameters are varied, include Hopf bifurcation, period doubling cascades 
leading to chaotic pulsations, Shilnikov-Hopf bifurcation etc. 

The new feature that was added to the setting is a variable initial concentration of the 
so-called deficient component which controls the reaction rate in the case of the combustion 
synthesis. We assume the concentration to be a periodically perturbed constant. While the 
original objective of our experiments was to study the dependence of the mean propagation 
velocity on the perturbation frequency and amplitude, we have noticed that the dynamics 
was qualitatively substantially different from that of the unperturbed problem. What we 
observed was a variety of quasi-periodic and complex periodic regimes with periods that 
are various multiples of the perturbation period. These observations present a convincing 
evidence that we are dealing here with the phenomenon of frequency locking. 

The physical phenomena modeled by periodically driven dynamical systems appear in 
many fields (see e.g. [TU]) such as lasers, superconductors (Josephson's junctions), mechanical 
engineering, etc. Theoretical as well as numerical studies show that periodic forcing can 
drive these systems to exhibit rich patterns of behavior that includes mode loking, with 
the mode-locked bands usually having the structure of the so-called Arnold's tongues in the 
amplitude-frequency parameter space. 

It is necessary to mention that the analysis or numerical simulation or a combination 
of thereof is usually carried out for a finite-dimensional system even if the original physical 
system is infinite-dimensional. The reduction in such cases is based on an assumption of a 
certain ansatz which consequently leads to an ODE. At the same time the perturbation is 
usually external to the base system and it is periodic in time. In our case, however, we make 
no attempt of such reduction, neither are we aware of any ansatz for a solution leading to it, 
while the internal forcing is periodic in space, and the system remains formally autonomous. 

Below we present some examples of numerical solutions and a crude map of a basic 
resonance band in the amplitude-frequency parameter space that represents an Arnold's 
tongue. We remind the reader that every point in the parameter space is a result of numerical 
simulation of a nontrivial free-interface problem for a partial differential equation. 

Thus, the sharp interface model of combustion synthesis presents a very natural and 
transparent example of the frequency locking phenomenon for a PDE. 

2 Free-interface problem 

In the context of combustion of condensed matter a wave of exothermic chemical reaction 
transforms a solid combustible mixture directly into solid product. The 1-D model involves 
differential equations for the temperature u of the mixture and the relative concentration of 
the so-called deficient reactant Z (see e.g. Shkadinsky et al. [llj): 

ut = nu^x + qW{Z,u), 
Zt = -W{Z,u), 

where k is the thermal diffusivity, W is the chemical reaction rate, and q is the heat release. 
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Due to a strong temperature dependence of the reaction there is a well defined narrow 
region (flame front) where the bulk of chemical reaction and the heat release occur. Thus the 
distributed chemical reaction can be replaced by the (5-function (see Zeldovich et al. [T2]). 



W = Z{s{t))g{u)6{x- s{t)) 



located at the interface x = s{t) between the fresh, Z = Z{s(t)) and burned Z = material. 

The equation with the 5-function source is rewritten as a system of two heat equations 
coupled at the interface. In the context of solidification with undercooling or rapid 
crystallization of thin films [S], the free-interface model studied below is conceptually even 
simpler: the latent heat of the phase transition released at the interface must be diffused by 
the surrounding matter. 

Therefore we shall be concerned with the following appropriately non-dimensionalized 
free-interface problem: find s{t) and u{x,t) such that 



and the derivatives and u~ are taken from right side and left side of the free interface 
respectively. At —oo the surrounding matter is assumed to be at the temperature of the fresh 
combustible, while the temperature of the burned matter as well as its gradient is bounded: 



Under reasonable assumptions on the kinetics functions g, one can rigorously prove that 
the free-interface problem ()2.1|) - ()2.4|) possesses global in time, uniformly bounded classical 
solutions (see j2]). One should not regard the rigorous proof of the existence and uniform 
boundedness of solutions as a futile academic exercise but rather view it as a verification of 
correctness of the model. It provides a firm foundation for the numerical simulations pre- 
sented in the paper. In particular, the proof of uniform boundedness of solutions underscores 
the dynamical robustness of the model. 

3 Overview of dynamics for unforced problem 

For reader's convenience in this Section we give a brief overview of the unforced dynamics 
(see jT]), i.e., with Zq{x) = 1. It is convenient to rewrite the non-equilibrium interface 
condition (|2.3p in the form: 



Ut = u^x, X ^ s{t), (2.1) 

u{x,0) =uo{x) >0, (2.2) 

Zoisit))g[uisit),t)]=v{t) fort>0, (2.3) 

u^{s{t),t) ~u^{s{t),t) =v{t) fort>0, (2.4) 




u{—oo,t) = 0, u{oo,t) < C, \ux{oo,t)\ < C. 



(2.5) 



V = g[u{s{t),t)] := 1 + aJ{u{s{t),t)) 



(3.6) 
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We shall assume that the function J{$,) = {g{^) — l)/a; is normalized in such a way that 



J(l) = 0, J'(l) = -1, (3.7) 

which can be achieved by rescaling variables. This assumption makes the problems with 
different kinetics identical in terms of linearization about the basic solution. We note that 
the variables are selected so that the interface propagates to the left. 

In order to clarify the meaning of the (positive) parameter a which is the main instability 
parameter, we note that for the Arrhenius type kinetics 

/ N r — 1) , , , 

" = «'"' = -°^P' .+ (l-^)». ' 

a is the scaled activation energy for the exothermic chemical reaction that occurs at the 
interface, and a is the temperature ratio of the fresh and burned material for the traveling 
wave solution. 

The problem ()2.1|) - ()2.4|) has a unique traveling wave solution 

exp{x + t), x<-t 
u, = { ^ ^ ^ _^ , = -t. (3.9) 

provided J is monotone. The linear stability analysis indicates that the loss of stability 
occurs via a supercritical Hopf bifurcation at Ocr = + 2, the corresponding frequency is 

LJcr — 1.03, Tcr = l-njuJcr — 6.1. 

We have demonstrated in P that depending on the parameters, the system exhibits 
dynamical regimes that include a Hopf bifurcation followed by classical Feigenbaum's cascade 
of period doubling leading to chaotic pulsations, a Shilnikov-Hopf bifurcation etc. 

For numerical simulations below we employed an explicit finite-difference numerical code 
in the coordinate system attached to the interface r] = x — s{t). Problem (j2.H) - (j2.4|) is 
considered on a finite interval —L < r] < L with the Dirichlet condition u{—L,t) = 

do 

simulating the decay of the solution at — oo and Neumann condition —{L,t) = reflecting 

orj 

stabilization of the temperature far behind the interface in the burned matter (solid phase) 
.In view of the fact that the dynamics of the problem requires a very fine temporal resolution 
(to resolve sharp changes in the velocity), our experience shows that there is no advantage 
in using implicit methods in this case (see [Tj for detail). 

We note that in mathematical terms, the free boundary problem in ()2.H) - ()2.4|1 governs 
the temporal evolution in the infinitely-dimensional "phase" space of functions u{x) and 
scalars v. One way to present results would be through graphs of u{x,t) versus x,t and 
time histories v{t) vs. t. We also represent dynamics through projections of the infinitely- 
dimensional phase space onto the 3-dimensional space [u{s(t) — 1, t) , u{s(t) ,t) , u{s(t) + 1, t)]. 
I.e., the functional profile u{.,t) is represented by three values: at the interface and at two 
points equidistant from it. 
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Table 1: Arnold's tongue of resonances (v) in the (a;, a)-plane; 0.5 < a; < 1.5 (horizontally), 
< a < 0.1 (vertically). 



5 



4 Frequency locking 



The main objective of this paper is to demonstrate the phenomenon of frequency locking 
for the driven free-boundary problem ()2.H) - ()2.4|1 . To this end we choose the initial mass 
concentration to be in the form of the harmonically perturbed unity: 



Note that unlike typical forced problems, the perturbation here is a function of the dependent 
variable of the problem s{t). Thus, strictly speaking, the problem remains autonomous. 
Therefore the driving frequency u and the frequency of the periodic solution, as seen in the 
power spectra and temporal dynamics in the figures below, are not the same. Of course, it 
is possible to introduce a new time r = s{t), which would make the perturbation "time"- 
dependent and the problem non-autonomous. However, it introduces a stiff highly-variable 
coefficient 1/f at the time derivative of the heat equation that makes numerical solution 
substantially more challenging. 

Parameters in ()3.8|1 were selected to be a = 0.1, a = 4.5 which correspond to a regime 
with simple periodic oscillations just past the Hopf biurcation for the unperturbed problem. 
We vary the forcing amplitude a from a = to a = 0.1 with an increment of 0.005 and the 
frequency u from u = 0.5 to a; = 1.5 with the step Auj = 0.05. The results show typical 
structures of the Arnold's tongues in the (a, a;)-parameter plane. Note that for every pair of 
parameters the computation involves numerical solving a nontrivial free-interface problem, 
and the computational cost of an exhausting analysis of an extended domain in the parameter 
plane becomes prohibitive. 

Table 1 represents a crude outline of one of Arnold's tongues. We note that in addition 
to the main tongue one can also observe what seems to be smaller resonant tongues that 
correspond to higher order resonances. The resolution of the computations however, does 
not allow us to make a definitive judgement concerning additional tongues. 

Next we discuss several representative solutions of the forced problem for various pa- 
rameter sets. Each of the Figures below contains a velocity profile (top), its power spec- 
trum (bottom- left), and a 3-dimensional projection of the orbit into the space [u{s(t) — 
l,t),u{s{t),t),u{s{t) + l,t)], where u{.,t) is the temperature (bottom-right). 

Fig. 1 depicts a simple periodic orbit. Fig. 2 shows a trajectory whose period corresponds 
to twice the period of the forcing (note that the forcing has a spatial not temporal period). 
Figs. 3 and 4 depict periodic orbits of high multiplicity in ultraharmonic and subharmonic 
domains respectively. Fig. 5 demonstrates an example of what appears to be a chaotic 
trajectory, which can be incurred from the power spectrum. This might indicate a presence 
of a period-doubling cascade. Finally, Fig. 6 shows a typical quasi periodic orbit. We remark 
that quasi periodic solutions with more than two basic periods have also been observed. 




(4.10) 



6 




7 




Fig. 3. Ultraharmonic resonance of multiplicity 4 for a = .08, = .65 
I 
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Fig. 4. Subliarmonic resonance of large multiplicity for a = .08, = 2.5 
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Fig. 5. Cliaotic orbit for a = .04, uj = 1.4 
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Fig. 6. Typical quasi-periodic orbit, a = .1, = 4 
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5 Concluding remarks 



Above we have demonstrated that the free-boundary problem ()2.1|l - ()2.4j) exhibits some basic 
behavior quite analogous to the prototypical finite-dimensional frequency locking systems 
such as, for instance, the driven van der Pol equation. Since the sole objective of this letter 
is to demonstrate the presence of the phenomenon of mode locking for the free-interface 
problem we have deliberately chosen a very limited setting, and, restrained ourselves from 
discussing a multitude of other questions that arise in the context of this complex dynamical 
phenomenon. 

There are a number of complex features observed in dynamical response of the van der Pol 
system j3 E] and other finite-dimensional counterparts that one would attempt to verify for 
our system . One would naturally ask, for instance, whether the Arnold's tongues represent 
a dense set with its complement forming a Cantor set in the (a, ci;)-plane,. whether the 
Feigenbaum sequences occur within the tongues etc. It should be noted that the dynamics 
of the unforced free-boundary problem varies dramatically as the nonlinearity parameter a 
increases. In this letter we considered only the simplest possible case corresponding to stable 
limit cycle following Hopf bifurcation. It would be extremely interesting to investigate how 
the frequency locking response changes with increasing a. We hope to be able to pursue 
these issues in the near future. 
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